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ABSTRACT 

Improvements of the Gaussian quadrature in 
conjunction with the Newton-Raphson iteration technique (TM 000 789) 
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It has been shown in the RB-71-35 that Gaussian quadrature sup- 
plemented by a Newt on-Raph son iteration technique is an effective method 
of calculation of tb^ Bivariate Nonnal (Tetrachoric) r • Subsequent studies 
made it apparent that significant improvements could easily be incorporated 
into the calculation with relatively little increase in complexity, or cost, 
of the computational technique. These improvements are discussed in this 
supplement. 

1. Estimates of the standard deviates h and k 
The value of k may be written in the form 

Ik - 

* ~ !c ^ ^ 

where, in the notation of the original paper, ^ (k) is equivalent to .5 
minus the marginal percentage, . In the original model, Hastings’ approx~ 
Imation without any modification was used to estimate k, although a remark 
was made that the result could be improved by an iteration technique. 

It is apparent from this form of the integral that we are faced with 
precisely the same problem as in our evaluation of r , namely, we must 
compute (or estimate) a variable upper limit of a definite integral. There 
Is no reason, therefore, not to use the same algorithmic technique, i.e. , 

Gaussian quadrature and Newton-Raphson iteration, to improve Hastings' esti~ 
mates. This seems even more feasible when one realizes the necessary calculation 
Ingredients, the Gaussian quadrature coefficients, the related weights, and 
the iteration structure are already available for the evaluation of r . 



2. The Gaussian Quadrature. 



The 5-point quadrature used in the original ^tiidy was quite rapid 
and gave acceptable values except where the joint and marginal values were 
close. However, an increase of only 3 points to an 8-point quadrature 
resulted in the convergence of many values which previously had failed. 

This increased accuracy is impor t not only in the evaluation of the 
final r integra.l, but of equal and perhaps greater benefit In establishing 
more accurate values of the h and k parameters which make up the function. 
Consider the following table: 



Different Quadrature Effects on h and k Calculation 



Area 


h 

true 


Hastings ’ 
Estimate 
Unmodified 


5-Point Quadrature 


8-Point Quadrature 


Value 


Iterations 


Value 


Iterations 


,5 


0 


-1.01- lO"^ 


-13 

-.3*10 


1 


-13 

-.4*10 


1 


.158655254 


1 


.999968 


1.0000004 


s 

2 


1.0000002 


2 


.022750132 


2 


2.000435 


2.000002 


2 


2.000001 


2 


.001349898 


3 


3.000314 


3.00022 


2 


2.999990 


2 



ThuS) a substantial improvement in the values of h and k is achie'^^ed with 
only two iterations, 

3m %e Starting Estimate 

In accordance with the above improvements, two terms of the series 
expansion were used instead of one, and the resulting quadratic equation im 
r solved to provide a better starting estimate. Extreme values again caused 
this estimate to exceed 1, consequently it was necessary to set limiting 
values as was done previously. Again, no one value seemed to assure convergence 
over the entire range of r . For example, a P value of .001131 (b =" 2, k = 3, 
r = .80) failed to converge with a starting estimate of .97 but converged to 




.8003 readily with a lower value. On the other hand, a P value of .477473 
(h = 0, k = 0, r = .99) failed with a starting estimate of .90 but con- 
verged in 5 iterations to .99096 with a starting estimate of .97. On the 
assimiption that the majority of r calculations will be within the range 
- .80 < r < .80 and only occasionally near the extreme values which tend 
to give the most computational difficulty, the bounds were set at + .80 
with a final pass using + .97 if the first fails to converge. Fairly 
extensive testing has resulted in the convergence of all "reasonable 
values by this method. 

4. The Convergence Criteria 

Two convergence values were used in the attached examples; 

1*10 for both h and k calculations and 1*10 for the r calculation. 

The effect of these values is, of course, evident in the above table. 

Summary 

Three versions I; ^ rithm thus for uset 

1) 5-point quadrature, unmodified Hastings' estimates of h and k. 

2) 5~point quadrature, improved estimates of h and k. 

J) E-^point quadrature, improved estimates of h and k. 

Ihe sttached sheet of computer output indicates the range oi; values for 
the 8-point quadrature. An average calculation for the 8-poitrt quadrature 
required .027 seconds per computed value of r vs ^ .020 foir n 5— point 

quadrecure. On 50 x 50 matrix of "live" data an a’/erage r rs-quired .018 
seconds using 8-point quadrature compared v;ith .013 seconds fo r 5 points. 

The stability of the ’dgher quadrature seems to justify its lEse. 
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